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Abstract — This paper studies the problem of determining the 
sensor locations in a large sensor network using relative distance 
(range) measurements only. Our work follows from a seminal 
paper by Khan et al. |1| where a distributed algorithm, known 
as DILOC, for sensor localization is given using the barycentric 
coordinate. A main limitation of the DILOC algorithm is that all 
sensor nodes must be inside the convex hull of the anchor nodes. 
In this paper, we consider a general sensor network without the 
convex hull assumption, which incurs challenges in determining 
the sign pattern of the barycentric coordinate. A criterion is 
developed to address this issue based on available distance 
measurements. Also, a new distributed algorithm is proposed 
to guarantee the asymptotic localization of all locaUzable sensor 
nodes. 

L Introduction 

Location based service is the fundamental issue in the area 
of sensor networks and it requires to solve the localization 
problem. The localization problem consists of two parts, 
namely, acquiring measurements and transforming them to 
coordinate information. In this paper, we consider the situation 
of using relative distance measurements only and focus on the 
localization procedure. 

Existing work on localization can be divided into two 
classes [2|: sequential methods and concurrent methods. A 
sequential method begins with a set of anchor nodes and 
computes the locations of other nodes one by one or group 
by group. A prominent example is the so-called trilateration 
method. Its advantage is that it is easy to implement, but it re- 
quires each location-unknown node to have three connections 
(edges) with location-known nodes, which is a sufficient but 
not necessary condition for localizability. 

A concurrent method starts with some initial estimate for 
the coordinate of every sensor node. Each node then updates 
its coordinate in a distributed or cooperative manner using 
the relative distance measurements with its neighbors and the 
estimates of the neighbors' coordinates. The iterative process 
terminates when the estimates converge, hopefully to the true 
coordinates. 

A novel concurrent method called distributed iterative 
localization (DILOC) was given by Khan et. al. in [l] 
based on the barycentric coordinate representation for sensor 



localization. The unique feature of this method is that the 
sensor locations can be expressed as a linear system, which 
can be solved iterative ly due to the desirable eigenvalue 
distribution of the linear system. Their method relies on two 
key assumptions on the network topology. First, all location- 
unknown nodes must be inside the convex hull of the anchor 
nodes. Second, every node, other than the anchor nodes, must 
be inside a triangle of three neighboring nodes. The main 
reasons of requiring these two assumptions are to avoid the 
difficulties in determining the sign patterns of the barycentric 
coordinate and to ensure the desired eigenvalue distribution of 
the hnear system so that the iterative algorithm asymptotically 
converges. These two assumptions are kind of restrictive, 
especially for a large sensor network when the sensing range 
for each node is limited and when the neighbors for each node 
cannot be arbitrarily arranged. 

In this paper, we follow the work of [ 1 1 by generalizing 
their method and eliminating the two key assumptions as 
mentioned above. The main idea is to employ a general 
form of barycentric coordinate representation which allows 
the coordinate of each node to be expressed as a linear 
function of the coordinates of any three neighbors. A criterion 
is developed, according to which the linear function can be 
determined using the relative distance measurements alone. 
The implication of this result is that all the sensor locations 
can be expressed as a linear system, just like in the standard 
DILOC case. However, the new linear system may not have 
a desired eigenvalue distribution like the standard DILOC 
algorithm to work. We then provide a new distributed iterative 
algorithm for localization. This is done by applying a diagonal 
pre-conditioner to the linear system. A simulation result is 
also provided to validate the effectiveness of our proposed 
algorithm. 

Notations: C denotes the set of complex number E denotes 
the Euclidean space. 1„ represents the n-dimensional vector 
of ones and /„ denotes the identity matrix of order n. The 
bold font of letter indicates vector and capital letter indicates 
matrix. A^jTj denotes a triangle formed by node i, j and k. 



II. Preliminaries and Problem Formulation 

A. Barycentric coordinates 

The barycentric coordinate, which was firstly introduced by 
August Ferdinand Mobius in 1827 [3J, is a geometric notion 
characterizing the relative position of one node with respect 
to its several neighbor nodes. For one node, say I with its 
Euclidean coordinate pi, and its three neighbor nodes, say i, 
j and k with their Euclidean coordinates pi, pj and pk in the 
plane, node Ts barycentric coordinate with respect to i, j and 
k is {au,aij,aik} satisfying 



Pi = aiiPi + aijPj + aikPk- 



(1) 



Especially, when aa + aij +aik — 1, the barycentric coordinate 
is called the areal coordinate because it can be expressed as 
a ratio of signed areas between specified triangles. As shown 
in Fig. [H the barycentric coordinate {an, aij, aik} is given by 



(2) 



where SAijk, SAiki, S^uj and SaijU are the signed areas 
of the corresponding triangles Aljk, Alki, Alij and Aijk. 
These areas can be calculated with pairwise internode distance 
measurements through Cayley-Menger determinant |,4J. For 
instance. 
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where dij, dik and djk are the distance measurements among 
node I, j and k, respectively. The sign of S^iik is positive if 
node / is on the left-hand side when one moves from node j 
to k, and negative otherwise. 




Fig. 1 . An illustrative example for the barycentric coordinate. 

To avoid the case that S/^ijk — 0, we need an assumption 
on the configuration of node /'s three neighbors as below. 
AO: For each node in the network, its three neighbors are not 
collinear. 

Note that, the computation of the coefficients in terms of 
dU, i.e., an, aij and aik, depends on the signed value of 



SAijk, SAiki and Sauj- If we only know the pairwise distance 
measurements, we can compute the square values of these 
areas and thus the absolute values of an, aij and aik- However, 
we cannot determine the signs of these areas and thus are not 
able to have the barycentric coordinate. 

B. Problem formulation 

The common used trilateration scheme for computing the 
coordinate of node I is to solve a group of equations like 



du = \\pi ~Pt\\ 
dij = lb; -PjW 
dik = \\pi -Pk\ 



(4) 



Here, p^, u G k, I}, is the Euclidean coordinate of node 
u and duv, u,v G {i,j,k,l}, is the distance measurement 
between node u and v. These equations can be solved in 
a sequential way if each node has at least three distance 
measurements to other nodes that know their coordinates. We 
call those sensor nodes, who initially know their coordinates, 
the anchor nodes. In contrast, we call the nodes, who do not 
know their coordinates initially, the normal sensor node or 
just sensor nodes for short. In this paper, a network under 
consideration is assumed to contain at least three anchor nodes, 
which is a necessity for uniquely localizing the network. 

Instead of solving these nonlinear equations in a sequential 
way. Khan et. al. provide an iterative algorithm |IT], named 
DILOC, to compute the coordinates of a network Q based 
on the barycentric coordinate presentation. 

Given a network Q, containing n nodes, its Euclidean 
coordinates can be written in a form like 



p = Ap, 



(5) 



where p £ C" is the aggregated Euclidean coordinate of Q 
and the nonzero inputs in the i-th row, i G {1, 2, • • • , n} are 
the barycentric coordinate of node i. For example, if node I 
has three neighbors i, j and k, then the l-\h row of A has 
nonzero entries in the positions corresponding to the i-th, j-th 
and fc-th columns. The other entries in the l-th row are all 
zeros. 

If we consider the first three rows of p to represent the 
positions of the three anchor nodes in Q, then in terms of the 
natural partition of anchor nodes and normal sensor nodes, we 
can partition A and p as 
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(6) 



where p^ and Ps correspond to the aggregate positions of 
anchor nodes and other normal sensor nodes, respectively. 
The nonzero entries of the l-\h row can be recognized as the 
weights of its three neighbors. So, all diagonal inputs of C are 
zeros. The matrix A is often treated as the adjacency matrix 
of the network. Thus, the representation of Ps can be written 
as 

p, = Cps + Bpa (7) 



and 



(/ - C)ps = Bpa. 



(8) 



According to eq. (Q, we could obtain an iterative algorithm 
to solve Ps under certain conditions. 



2,{t + I) = Cz,{t) + Bpait) 



(9) 



where Zg G C"^'^ represents the estimate of the coordinate 
of the normal sensor nodes in Q. 

In H], a sensor network to be localized is assumed to satisfy 
the following two assumptions. 

Al: All sensor nodes lie inside a convex hull formed by the 
anchor nodes. 

A2: Each sensor node / lies inside a convex hull formed by 
its three neighbors. 
In the 2D case, a convex hull is actually a triangle formed 
by three nodes. Thus, assumption A2 leads to two constraints 
of the barycentric coordinate, i.e.. 



an + aij + aik = 1, 
< au,aij,aik < 1. 



(10) 
(11) 



In m. Khan et. al. proved that the spectral radius of C is 
less than 1 when (fTOl i and (fTTT i hold. Thus, the estimate of 
the coordinate in (|9]l can converge to the true value in the 
Euclidean coordinate system. 

In this paper, we will relax the assumptions Al and A2 and 
address a distributed algorithm to compute the locations of 
sensor nodes no matter whether they lie inside a convex hull 
or not. After relaxing these two assumptions, two problems 
need to be addressed. First, how to determine the signs of 
the barycentric coordinate when one node lies outside the 
convex hull of its neighbors. Second, when the convex hull 
assumption is dropped, the matrix C in system ^ might not be 
Schur Then, how to provide a convergent iterative algorithm 
to compute the coordinate. 

III. Sign pattern determination for the 

BARYCENTRIC COORDINATE 

For u,v G {i,j,k,l}, we use (t„„ G {1,-1} to indicate the 
sign of Quv Suppose node / is locahzable. It is known that 
no matter a node I lies inside the convex hull of its three 
neighbors or not, the barycentric coordinate {aii,aij,aik} 
obtained from eq. (|2|i must satisfy an + aij + ait — 1. Thus, 
given \aii\, |a;j|, \aik\, which can be calculated from (|3]l, the 
problem of determining the sign pattern of the barycentric 
coordinate is equivalent to solve the following equation 



<7ii\aii\ 



(yij\aij\ 



crik\aik\ 



(12) 



where (Jii,<Jij and <Tik take values either 1 or —1. 

In the following we will discuss whether (fT2] | has a unique 
solution. Moreover, if au^aij and aik can not be uniquely 
solved from ( fT2] i. we then explore other range based conditions 
to determine the sign pattern. 

The first case that (fT2] l does not have a unique solution is 
that one of \aii\, \aij\, \aik\ equals to zero. That is, a node I 



lies on the line aligned with one of three edges of the triangle 
formed by its three neighbors, according to (|2). Without loss 
of generality, say an — 0. For this case, an can be either 1 
or —1. But the other two signs aij and aik can be determined 
according to the following criterion. 

{{an, 1, 1} if \aij\, \aik\ < 1, 
{an, 1, -1} if \aij\ > 1, \aij\ > \aik\, 
{an, -1, 1} if \aik\ > 1, \aik\ > \aij\. 

(13) 

If node I does not lie on the boundary lines, there are totally 
7 possible sign patterns as the pattern { — 1,— 1,— 1} is not 
possible due to (fT2l i. According to ^ and the definition of the 
signed areas, the seven possible sign patterns of {an, aij, aik} 
are shown in Fig. |2] In the following lemma, we characterize 



{1,-1,1} 




1-1,-1,1} 



Fig. 2. Seven possible sign patterns for {(^n, ryij , crik}. 

the second case when the sign pattern can not be uniquely 
solved from (fT2l i. 

Lemma 1: Given \an\ ^ 0, \aij\ ^ 0, and \aik\ ^ 0, 
the solution of (fT2l l does not result in a unique sign pattern 
{ch, Cij, aik{ if and only if one of them, saying an, satisfies 
\an\ = 1, and |a;j| = \aik\. 

Proof: (Sufficiency) If \an\ — 1 and \aij\ = \aik\, it 
can be inferred from (IT2t that {an, aij, era} = {1,1,-1} or 
{ch, fij, aik{ = {1, —1, 1}. That is, (fT2] i does not result in a 
unique sign pattern. 

(Necessity) Suppose there are two sign patterns both satis- 
fying (fTST i. That is, it holds that 



[\an 



\aik 



[1 1] 



(14) 



where vi,V2,V3 G {[1 1], [-1 1], [-1 - 1], [1 - 1]}. This 
means a positive combination of vi, V2, V3 equals to [1 1] (see 
Fig. [3]i. Consequently, there must be [1 1] for one of Vi, V2, 
and V3. Without loss of generality, we assume vi = [1 1]. 
Next, we consider different choice of V2. If V2 equals to [1, 1] 
or [—1,-1], we will have V3 equal to [—1 — 1] or [1 1]. In 
this way, two solutions {an, aij, aik} are actually identical. If 
V2 equals to [—1 1] or [1 — 1], V3 must equal to [1 — 1] or 
[—1 1]. Then according to ( fT2] i. we know that \an\ ~ 1 and 
\aij\ = \aik\. ■ 
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Fig. 3. An illustration for the necessity proof. 



Next, we present a result on how to determine the sign 
pattern using the range based information when it can not be 
uniquely solved from (fTSI i. 

Lemma 2: Given \aii\ = 1 and \aij \ = \aik\ ^ 0, suppose 
Zijk is an acute angl^H 

1) CTji = — 1 if and only if 

dji = dik, dki = dij, and d| = 2d^j + 2d,\. - d^j,. 

2) If an = -1, then {(7;^, ctjj , era} = {-1,1,1}. 

3) If (Tii = 1 and < df^ + d|, then {(Jii,aij,aik} = 
{1,1,-1}. 

4) If an = 1 and > df^ + d^, then {(Jii,aij,aik} = 
{1,-1,1}. 

Proof: 1) (Necessity) If an = —1, then we have an = 
— 1. Moreover, since \aij\ — \aik\, it follows from (fT2] i that 
aij = aik ~ 1. Thus, SAijk — ~SAijk, Sauj = SAijk and 
SAiki — SAijk- Comparing the sign pattern with the ones 
described in Fig. |2] we know that the only option of I is at the 
location of / in Fig. |4l which forms a parallelogram together 
with nodes i, j and k. Hence, we can obtain directly that dji ~ 
dik, dki = dij. Furthermore, according to the parallelogram 
law, we have 

4=24+2^2^-4,,. (15) 

(Sufficiency) If dji — dik and dki = dij, we can draw two 
circles centered at j and k with radius dik and dij, respectively. 
These two circles will have two intersection points. One of the 
two intersection points is I and we denote the other by /*. 
From the necessity proof, it is known that when node I is at 
the location of I , it satisfies (fTsT l. On the other hand, we will 
show that when node I is at the location of /*, it does not 
satisfy (flSl l. (To see this, it remains to show that d^i'" ^ dn*. 
Suppose by contradiction that d^i"i — dn-. Then, recalling the 
fact dj^ni = dji* and d,;'" = dki*, we have that nodes i, j 
and k are on the perpendicular bisector of the line segment 
I I* and so they are colinear, a contradiction to assumption 
AO.) Therefore, it can be concluded that due to the condition 
d^i = 2d^j + 2d^f. — 4j,, node I must lie at the location of I . 
Thus, according to the sign patterns described in Fig. |2] we 
can obtain that an = —1. 

'if Zijk is not acute, then Zikj must be acute and the conditions in the 
lemma can be modified accordingly. 



l' i l" 




V" 

Fig. 4. An example of Aijk and node I. 



2) If an = -1, then {an,aij,aik} = {-1,1,1}, which is 
shown in the necessity proof of 1). 

J) If an — 1, then we have an — 1. For this case, SAijk — 
SAijk according to (|2]i. Therefore, node I must be on the line 
that is parallel to the edge jk and crosses node i. On this line, 
there are two nodes, saying node I and / as shown in Fig. |4] 
whose distances to node i are equal to dn. 

For the triangle /S.ijl , according to the cosine law, it holds 
that 

d^.^, = 4 + 4' ^ '^'^i'j'^ii' cos Zjil' . 

Since I I is parallel with jk, we have Zjil ~ Zijk. Then 

we know (P:,, < 4 + d^,, because Zijk is an acute angle. 

]i V ii ° 

Similarly, for the triangle Aijl , we could obtain that d^^,, > 

4 + 4„. 

Therefore, if 4; < dfj + d\, then node / must be at the 
locatio of Thus, according to the sign patterns described in 
Fig. ID we obtain that {an, <Jij,mk} = {1, 1, -1}. 

4) Following the argument in 3), we know that if 4; > 
dfj + 4;, then node / must lie at the location of / . Then 
again from the sign patterns described in Fig. |2] we obtain 

that {(7;i,CT;j,crifc} = {1,-1,1}. ■ 

Finally, we summarize the above results to provide an 
algorithm of determining the sign pattern based on the range 
measurement information. The pseudo code is given in Algo- 
rithm [T] 

IV. Distributed localization algorithm 

In this section, we develop a convergent iterative algorithm 
for localization, that is, to solve from (|7]i. Notice that for a 
general barycentric coordinate, the spectral radius of C may be 
larger than 1 . It means that the iterative form presented in ^ 
may not converge. In the sequel, we will modify it to have a 
new convergent iterative algorithm to compute the coordinates 
of localizable sensor nodes. 

Pre-multiplying a diagonal matrix K on both sides of (|7]l, 
we can obtain 

Kp, = KCps + KBpa. 



Algorithm 1 Determining the sign pattern of node Vs barycen- 

tric coordinate. 

Input: \aii\, \aij\, \aik\, du, dij, dik, dij, dik, djk- 
Output: {(Tii,(Tij,aik}- 
Solve eq. ([T2]i 

if the solution is unique then 

Return au,aij,aik 
else if one of |a/i|, \aij\, \aik\ equals to then 

Determine {(7;^, cr;^ , crifc} according to (flJl l 
else if dji = dik, du = dij and dfi = 2dfj + 2djf, - d^^, 
then 

{an, (J I j, (J Ik} = {-1,1,1} 
else if d^, < + d| then 

{crii,<7ij,aik} = {1,1,-1} 
else if d^, > df^ + d| then 

{crii,<7ij,aik} = {1,-1,1} 
end if 



Adding Ps on both sides of the above equation and packing 
the terms, we have 

p, = iI-K{I-C))ps+KBpa. 

Recall that Pa is constant. So we could consider the following 
iterative algorithm 

z,{t + 1) = (/ - K{I - C))z,(<) + KBpa, (16) 

where is the estimate of the coordinates of the normal 
sensor nodes. Here, the key is to find an appropriate diagonal 
matrix K so that / — K{I — C) is Schur. 

The above iterative form of localization can be implemented 
in a distributed way, that is, 

z,{t + I) = z,{t) ~ kA z,{t) - T.a,jZj{t)\ (17) 

where Zi is the estimate of sensor node i's coordinate, ki is 
the ith diagonal entry of K, and a^j, j G A/i, is the barycentric 
coordinate of node i with respect to its three neighbors. 

However, finding a diagonal pre-conditioner K in a com- 
pletely distributed way to make I — K{I — C) Schur is a 
challenging task. Next we introduce a cluster-based approach 
for the design of K, which is partially distributed and does 
not require to collect all the information of the whole network. 
For a sensor network Q, if the anchor nodes are not far away 
from each other, then suppose Q can be partitioned into a 
set of clusters Cyo, Qi, • • • , Gm such that Qq is the cluster 
of anchor nodes, and for any cluster Qs, s £ {1, • • • , ?t7,}, 
the neighbors used to define the barycentric coordinate of 
any node in Qs belong to t/o U • • • U C/s- This is the case 
when the set of clusters is sequentially localizable. It is 
common in practice as the neighboring topology is usually 
dependent on the configuration in the Euclidean space and 
limited communication ranges. How to detect such a set of 
sequential clusters is of independent interest. Readers may 
refer to [8] and g]. 



For each cluster Qs, s — 1, . . . , m, we propose a scheme 
to design Kg (the corresponding block of K for the cluster 
Qs). The idea is inspired by the work of [5] and [6 | using the 
continuity property of eigenvalues with respect to the values 
of Ks- Denote by Ls the square block submatrix in / — C 
corresponding to the cluster Qs (taking the rows and columns 
of I—C indexed by the nodes in Qs). Suppose the total number 
of nodes in Qs is n^. We denote Ks — diag(fcsi , fcs2 1 ■ • • i fcs„ ) 
and denote by L^, j — 1, . . . , Us, the jth principal submatrix 
of Ls- The procedure of designing Ks for each cluster Qs is 
summarized in Algorithm |2] 



Algorithm 2 Finding the diagonal pre-conditioner Ks for 
cluster Qs- 
1: Input: Ls 

2: Output: Ks = diag(fcsi , fc^a , • ■ • , fcs„, ) 
3: for j — I : Us do 

4: Find k'^. such that the eigenvalues of 
diag{fc^^ ,kg. }Li lie in the open right-half- 
plane. 

5: end for 

6: Find sufficiently small e > to make the eigenvalues of 
ediag{fcg , . . . , fc^ }Ls inside the unit circle centered at 
(1,0). 

7: return fc,^. = efc^^, j = 1, . . . , n^. 



Next, we discuss why Algorithm |2] can succeed in finding 
the diagonal pre-conditioner Ks- According to the loop in the 
algorithm, it is certain that an appropriate k'^^ can be found 
first so that k'^^L\ is in the open right-half-plane. Thus, the 
eigenvalues of 



are k'g_^L\ and 0. By the continuity property of eigenvalues, 
we can then find k'^^ in the neighborhood of the origin such 
that the eigenvalues of 

both lie in the open right-half-plane. Repeating the argument 
leads to the finding of a set of k'^., j = 1, . . . , n^, such that 
the eigenvalues of diag{/c^^ , • • • , K„ all lie in the open 
right-half-plane, for which, a small e > can then be chosen 
to shrink the eigenvalues of ediag{fc^^ , • ■ • , K„ inside the 
unit disk centered at (1,0). 

To implement the algorithm of finding a diagonal pre- 
conditioner for the cluster Qs, a randomly selected node in 
the cluster acts as a cluster head and collects the barycentric 
coordinates of nodes in the same cluster, i.e., Lg- It computes 
an appropriate Ks according to Algorithm |2] and then sends 
ks to the individual nodes in the cluster. If the cluster size is 

Sj 

medium, the required communication cost is acceptable. 

After obtaining the diagonal pre-conditioners for all clusters, 
the localization algorithm (fTSI l is fully distributed, requir- 
ing only the exchange of the estimate from the neighbors. 



Moreover, it can be known that the iterative algorithm (ITS] 
is globally asymptotically convergent due to its linear form, 
while most existing localization work (e.g., |7|) based the sub- 
gradient method only ensures local convergence. 

V. Simulation 

In this section, a sensor network with 12 nodes is consid- 
ered. As shown in Fig. |5(a)[ three anchor nodes are connected 
by black lines and other nine sensor nodes are marked by red 
stars. The blue lines with arrows represent the neighboring 
topology in localization. In this example, both assumptions 
Al and A2 are not satisfied. 

The 12 nodes are grouped into four clusters, i.e., Qq ~ 
{1,2,3}, Gi = {4,5,6}, 02 = {7,8,9}, and 03 = 
{10,11,12}. The barycentric coordinate of each node is 
calculated according to Algorithm [T] based on the range 
measurement information. The diagonal pre-conditioner Kg, 
s = 1,2,3, is obtained utilizing Algorithm |2] The coordinate 
of each node is then iteratively calculated in terms of (fTTl i in a 
distributed way. The trajectories of the estimates are shown in 
Fig. |5(b)[ from which it is seen that the estimates asymptoti- 
cally converge to the true coordinates from an arbitrary initial 
guess. 

The residual error defined plotted in 

Fig. |5(c)| with respect to the iteration steps. Though the 
localization algorithm is executed in parallel, the convergence 
process takes three stages, as being observed in Fig. |5(c)| be- 
cause the localization of a cluster depends on the localization 
of the cluster closer to the anchor nodes. 

VI. Conclusion 

In this paper, we develop a distributed algorithm to com- 
pute the locations of sensor nodes based on the barycentric 
coordinates. Two critical problems are solved, leading to the 
success of globally convergent localization. First, for a general 
configuration that does not require every node to be inside 
a convex hull of its neighbors, a range information based 
algorithm is proposed to determine the signs of the barycentric 
coordinates and therefore the barycentric coordinates. Second, 
a distributed iterative algorithm is obtained with the global 
convergence ensured diagonal pre-conditioner designed based 
on a partially distributed cluster scheme. Future work includes 
analysis of localizability and convergence rate of the proposed 
approach, and localization performance in the presence of 
measurement noises. 
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